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Following recent claims relative to the question of large anisotropy production in regular bouncing 
scenarios, we study the evolution of such anisotropies in a model where an Ekpyrotic phase of 
contraction is followed by domination of a Galileon-type Lagrangian which generates a non-singular 
bounce. We show that the anisotropies decrease during the phase of Ekpyrotic contraction (as 
expected) and that they can be constrained to remain small during the non-singular bounce phase (a 
non-trivial result). Specifically, we derive the e-folding number of the phase of Ekpyrotic contraction 
which leads to a present-day anisotropy in agreement with current observational bounds. 
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I. INTRODUCTION 



The idea of a non-singular bouncing universe as an al- 
ternative to the singular Big Bang paradigm has been 
attractive on philosophical grounds for a long time. 
More recently [T] [2] it has been realized that a non- 
singular bouncing cosmology with a matter-dominated 
initial phase of contraction - a "matter bounce" cosmology 
- may provide an alternative to the inflationary paradigm 
of cosmological structure formation (see e.g. [3] for re- 
cent reviews). Since the curvature perturbation £ grows 
on super-Hubble scales in the contracting phase, fluctu- 
ations on long wavelengths get boosted by a larger fac- 
tor than those on small scales. In a matter-dominated 
phase of contraction the relative enhancement factor is 
precisely such as to transform an initial vacuum spectrum 
into a scale-invariant one. Studies in many non-singular 
models |3] have shown that the scale-invariance of the 
spectrum of C 1 before the bounce is maintained through 
the bounce on length scales which are larger than the 
duration of the bouncing phase (the phase in which the 
violation of the usual energy conditions is localized) 2 . 
Hence, the matter bounce scenario leads to the same 
shape of the spectrum of primordial curvature fluctua- 
tions on super-Hubble scales as inflationary cosmology. 
Starting with vacuum initial conditions, the fluctuations 
are also Gaussian and coherent. A specific prediction of 
the matter bounce with which the scenario can be dif- 



* Electronic address: yifucai@physics.mcgill.ca 
^Electronic address: rhb@physics.mcgill.ca 

* Electronic address: peterOiap. 

1 The variable £ is the curvature fluctuation in comoving gauge. 
It is proportional to the canonical fluctuation variable v - the 
Sasaki-Mukhanov variable 5,6 - in terms of which the action 
for cosmological fluctuations has canonical form. See e.g. [7] for 
a comprehensive review and [8] [9] for introductory overviews of 
the theory of cosmological perturbations. 

2 Note that there are counterexamples as discussed in I10| and 
more recently in 



ferentiated from inflation is the shape of the bispectrum 

Bouncing cosmologies, however, face several problems. 
First of all, new physics is required to provide the bounce 
(a pure GR bounce sourced by a regular fluid also leads 
to dynamical instabilities [13], while with a scalar field 
a positive spatial curvature [14] or non-canonical kinetic 
term |15| is required). Such new physics can either come 
from introducing new forms of matter such as phantom or 
quintom fields [IB 1 , ghost condensates [TTlfTB] . Galilcons 
|19j . S-branes [2D] or effective string theory actions |21| . 
It can also arise from corrections to Einstein gravity, as 
e.g. in the non-singular universe construction of |22| . 
specifically applied to the bouncing scenario in |23| , in 
the ghost-free higher derivative gravity model of [23], in 
terms of torsion gravity [25] , or in Hof ava-Lifshitz gravity 
|26| . Finally, a bouncing phase can also originate from 
pure quantum cosmological effects [57]; those also pro- 
vide a natural framework in which a dust-dominated con- 
traction is easily implemented to yield a scale-invariant 
spectrum of perturbation |28| . 

A more serious problem for bouncing cosmologies is 
the fact that the contracting phase is unstable to various 
effects. As discussed e.g. in [25], many matter bounce 
models are unstable to the addition of radiation. This 
arises since the energy density of radiation scales as a -4 , 
where a(t) is the cosmological scale factor, whereas the 
density of matter grows as a~ 3 . In some models (e.g. 
|16| ) this will lead to a "Big Crunch" singularity instead 
of a smooth bounce. Some constructions (e.g. those of 
[181 [22] , [2D] and [3D] ) are free from this problem. How- 
ever, a problem for all constructions mentioned so far is 
the instability to the growth of anisotropic stress, whose 
associated energy density grows as a -6 . This is the fa- 
mous BKL instability of a contracting universe |31j . 

A solution of this anisotropy problem can be realized in 
the Ekpyrotic scenario |32) . in which it is postulated that 
a matter field with equation of state parameter w 3> 1 
(where w = p/p is the ratio of pressure p to energy den- 
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sity p) is dominant in the contracting phase 3 . Such an 
equation of state can be realized by treating the dominant 
form of matter as a scalar field with negative exponen- 
tial potential. Since the energy density of the dominant 
matter then scales with aT q with g > 6, anisotropics be- 
come negligible and the BKL instability is avoided [37_ 4 . 
In a recent paper |38| . a subset of the present authors 
introduced a scalar field with an Ekpyrotic potential to 
construct a matter bounce scenario which is free from the 
BKL instability problem. 

The Ekpyrotic scenario in its original formulation [32] 
involves a singular bounce. In addition, the curvature 
spectrum of £ is an n s — 3 spectrum rather than a scale- 
invariant n s = 1 one [39-42J. Hence, without non-trivial 
matching of £ across the bounce, one cannot obtain a 
scale-invariant spectrum at late time 5 . To solve this 
problem, a new and non-singular version of the Ekpyrotic 
scenario [IS] was proposed in which a second scalar field 
is introduced which does not influence the background 
dynamics but develops a scale-invariant spectrum which 
starts out as an isocurvature mode but which is trans- 
ferred to the adiabatic mode during the evolution. The 
second field can also be given a "ghost condensate" La- 
grangian |47| in which case it mediates a non-singular 
bounce. However, as has been pointed out in [55], in 
this "New Ekpyrotic" scenario the anisotropics which are 
highly suppressed during the contracting phase again 
raise their head and lead to a BKL instability. 

In our previous work [38 , we argued qualitatively that 
in the model we considered the anisotropics remained 
negligibly small during the bouncing phase. The reason 
for the difference compared to what happens in the model 
of |46| is that in our model the kinetic condensate which 
grows as the bounce is approached does not need to de- 
crease again by the time of the bounce point. This leads 
to a shorter bounce time scale and to different dynamics. 

In this paper we carefully study the development of 
anisotropics in the bouncing cosmology with an Ekpy- 
rotic phase of contraction introduced in [38J. We work 
in the context of a homogeneous but anisotropic Bianchi 
cosmology in which the scale factors in each spatial di- 
mension evolve independently. We are able to show that 
no BKL type instability develops, in agreement with 
what the study of [38 indicated. Our work thus shows 
that the arguments against non-singular (as opposed to 
singular) bouncing cosmologies put forwards in |48| do 



3 There are other approaches to address the anisotropy prob- 
lem. For example, nonlinear matter terms may smooth out the 
anisotropies |33| . Adding quadratic R a pR a P terms to the grav- 
itational action can also prevent the BKL instability 1341 . 

4 Note, however, that including anisotropic pressures may reintro- 
duce instabilities towards anisotropy generation |35] . 

5 However, the spectrum of the Bardeen potential <& is scale- 
invariant [43 j . and, as argued in |10| and shown explicitly in some 
examples 44 45 , it is this spectrum which may pass through 
the bounce, thus yielding a scale-invariant spectrum of curva- 
ture fluctuations at late times. 



not apply to all non-singular bouncing cosmologies. 

The outline of this paper is as follows. In the next sec- 
tion we review the bounce model introduced in |38| and 
derive the resulting equations of motion for a homoge- 
neous but anisotropic universe. In Section 3 we analyt- 
ically study the background dynamics in each phase of 
the cosmological evolution from the initial matter phase 
of contraction through the Ekpyrotic phase to the bounc- 
ing phase and the subsequent fast-roll expanding period. 
Specifically, we determine the decay or growth rates of 
the anisotropy parameter in each phase. In Section 4 we 
solve the dynamical system numerically and present our 
final results. We close with a general discussion. 

A word on notation: We define the reduced Planck 
mass by M p] = 1/^8ttG n where G N is Newton's gravi- 
tational constant. The sign of the metric is taken to be 
(+, — , — , — ). Note that we take the value of the mean 
scale factor at the bounce point to be a B — 1 throughout 
the paper. 



II. A NONSINGULAR BOUNCE MODEL 

We consider a nonsingular bounce model in which the 
universe is filled with two matter components, a cosmic 
scalar field cf> and a generic matter fluid, as proposed in 
Ref. [55J (which, in turn, is based on the theory devel- 
oped in [49] ) . The Lagrangian of <fi is given by 

C[<f>(x)]=K(<i>,x) + G{<j>,x)n4>, (l) 

where K and G are functions of (p and its canonical ki- 
netic term 

X = l -d^d^, (2) 
while the other kinetic terms of </> include the operator 
□0 = gTVuVut- (3) 

Variation of the above scalar field Lagrangian mini- 
mally coupled to Einstein gravity leads to the following 
corresponding energy momentum tensor 

T+, = {-K + 2XG^ + G. x V (y XV a 4>)9^ 
+{K,x + G, A -D0 - 2G, )V^V„0 
-GxtV.lV^ + V.lV^), (4) 

in which we use the notation that JF^ and T x denote 
derivatives of whatever functional F(4>, X) may be with 
respect to <f> and X, respectively. 

For the model under consideration we choose: 

K{4>, X) = Ml [1 - g(<j>)} X + pX 2 - V{4>), (5) 

where we introduce a positive-definite parameter /3 so 
that the kinetic term is bounded from below at high en- 
ergy scales. Note that the first term of K involves 
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since in the present paper we adopt the convention that 
the scalar field <fi is dimensionless. 

The function g((f>) is chosen such that a phase of ghost 
condensation only occurs during a short time when <j> ap- 
proaches = 0. This requires the dimensionless function 
g to be smaller than unity when \<f>\ ^> 1 but larger than 
unity when approaches the origin. To obtain a nonsin- 
gular bounce, we must make an explicit choice of g as a 
function of <f>. We want g to be negligible when \4>\ 1. 
In order to obtain a violation of the Null Energy Con- 
dition after the termination of the Ekpyrotic contracting 
phase, g must become the dominant coefficient in the 
quadratic kinetic term when (f> approaches 0. Thus, we 
suggest its form to be 



9(4) 



2.9o 



(0) 




where go is a positive constant defined as the value of g 
at the moment when = 0, and is required to be larger 
than unity, go > 1. 

We have also introduced a non-trivial potential V for 0. 
This potential is chosen such that Ekpyrotic contraction 
is possible. It is well known that the homogeneous tra- 
jectory of a scalar field can be an attractor solution when 
its potential is an exponential function. One example is 
inflationary expansion of the universe in a positive- valued 
exponential potential, and the other one is the Ekpyrotic 
model in which the homogeneous field trajectory for a 
negative exponential potential is an attractor in a con- 
tracting universe. For a phase of Ekpyrotic contraction, 
we take the form of the potential to be 



F(0) 



2Vn 



(7) 



where Vq is a positive constant with dimension of (mass) 4 . 
Thus the potential is always negative and asymptotically 
approaches zero when \(f>\ 3> 1. Ignoring the second term 
of the denominator, this potential reduces to the form 
used in the Ekpyrotic scenario |32| . Both functions g(0) 
and V{(f) are shown on Fig. [ljwith the parameters used 
in the later parts of this work. 

The term G(cf>, X) is a Galileon type 6 operator which 
is consistent with the fact that the Lagrangian contains 
higher order derivative terms in 0, but the equation of 
motion remains a second order differential equation. Phe- 
nomenologically, there are few requirements on the ex- 
plicit form of G(<fi,X). We introduce this operator since 
we expect that it can be used to stabilize the gradi- 
ent term of cosmological perturbations, which requires 
that the sound speed parameter behaves smoothly and is 
positive-definite throughout most of the background evo- 
lution. For simplicity, we will choose G to be a simple 



See 1361 for a discussion of Galileon type Lagrangians. 



Figure 1: Model functions g((f>) and V(cf>) as given by Eqs. Q 
and J7J, with background parameters taken as for the follow- 
ing evolution figures, namely as in Eqs. (65 1 and (66 1. 



function of only X: 



G(X) = yX, 



(8) 



where 7 is a positive-definite number. 

We now turn to the study of the cosmology of this 
model. In order to characterize a homogeneous but 
anisotropic universe, we take the metric to be of the form 



d^ = dt 2 



V^^Te^VV 



(9) 



where t is cosmic time, a % are linearly independent at 
all points in space-time and form a three dimensional 
homogeneous space. 

In the case of a Ricci flat space, one can consider the 
projection o % = da; 1 and thus the metric is of Bianchi 
type-I form. The factor a(t) can be viewed as the mean 
scale factor of this universe, and the functions e® 1 ^ de- 
scribe the correction of anisotropics to the scale factor. 
Since the values of scale factors can be re-scaled arbitrar- 
ily, one can impose an additional constraint 



5> = o. 



(10) 



Then, one can immediately define a mean Hubble param- 
eter as follows, 



a 



(11) 



and the individual Hubble parameters along spatial di- 
rections are given by, 



H, 



1 d 

ae 0i dt 



[ae 6i ) = H + 6 L , (no sum) 



(12) 
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where the overdot denotes the derivative with respect to 
cosmic time t. 

Since we are interested in studying anisotropics rather 
than inhomogeneities we can treat the matter fields to be 
homogeneous, which implies (f> is only a function of cosmic 
time. Thus, the kinetic terms of the homogeneous scalar 
field background become 



X = 



□0 = 



1 ;> 



3H(f>, 



(13) 



so that, for this background, the energy density of the 
scalar field is 

P4> = \m 2 {1 - g)4? + + 3jH^ + V(cj>), (14) 



and the pressure is 



1 



M 2 (l-g)4> 2 + 



^ 2 4>-V(0), (15) 



as follows by computing the diagonal components of the 
stress-energy tensor Q. 

Additionally, the matter fluid contributes its own en- 
ergy density p m and pressure p m , and usually they are 
associated with a constant equation-of-state parameter 
w m = Pm/ Pm- Namely, for normal radiation, w m = |, 
while for normal matter, w m = 0. 

To derive the equation of motion for (f>, one can either 
vary the Lagrangian with respect to <f> or, equivalently, 
require that the covariant derivative of its stress-energy 
tensor vanishes. This yields 



V4> + V<j> + = 0, 
where we have introduced 



V 



V 



{l-g)M 2 + 6T ff0 + 3/?^ 2 + 



(16) 



(17) 



3(1 - g)M 2 H + \9yH A - ^M^g t4> ) + 30H<f? 



37 2 
2M 2 



3G 



x 



2M 2 



(pm + Pi, 



(18) 



From Eq. (16 1, it is clear that the function V determines 



the positivity of the kinetic term of the scalar field and 
thus can be used to determine whether the model con- 
tains a ghost or not at the perturbative level; the function 
T> on the other hand, represents an effective damping 
term. By keeping the first terms of the expressions of 
V and V and setting g — 0, one can recover the stan- 
dard Klein-Gordon equation in the FRW background. 
Neglecting the other terms is a good approximation when 
the velocity of <fi is sub-Planckian. Note that the friction 
term T> contains the contributions from anisotropic fac- 
tors and matter fluid, which can be suppressed for small 



values of <j>. However, these terms will become important 
during the bouncing phase where 4> reaches a maximal 
value. For simplicity, in the following we will consider 
matter fluid is cold and thus w m = 0. 

Finally, we can write down Einstein equations in this 
background, given by 



Mi 



"■pv 2 



_ rpc/> 



T, 



(19) 



Once expanded in components, this tensor equation 
yields the effective Friedmann equations, 



H 2 = 



H 



Pt 

3M|, ' 6 

Pt +Pt _ 

2M 2 2 ^ 



(20) 
(21) 



where p T and p T represent the total energy density and 
pressure in the Bianchi type-I universe, i.e., the sum of 
the contributions of the scalar field and the fluid. 

Moreover, combining the spatial component of Ein- 
stein equation with the constraint equation (10 1 yields 



6i + 3H9 l = 0, 



from which it follows that 



i(t) = M e ,i 



(22) 



(23) 



where a B is the mean scale factor of the universe at the 
bouncing point. The coefficients Mg^ are integral con- 
stants with a dimension of mass. According to the con- 
straint equation (fl0| , one can read off that 



(24) 



Plugging Eq. ( 23 1 into Eq. ( 20 1 shows that one can 



introduce an effective energy density of anisotropy 



M 2 

Fl 

2 



(25) 



whose evolution as 1/a 6 implies an effective equation-of- 
state parameter equal to Wg = 1. We see that this effec- 
tive energy density increases faster than that of pressure- 
less matter or radiation in a contraction universe. This 
is the source of the BKL instability of the contracting 
phase of many bouncing cosmologies. 



III. BACKGROUND EVOLUTION 

The initial conditions of our model are chosen (as in 
|38| 1 such that we start in a contracting phase dominated 
by regular matter. Since the energy density of the Ekpy- 
rotic scalar field <p grows faster than that of regular mat- 
ter, <fi will at some time begin to dominate the energy 
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density. At this point, the Ekpyrotic phase of contrac- 
tion begins, and lasts until the nonsingular bounce in- 
terval begins (this is the phase where the new physics 
effects dominate) , followed by a period of fast-roll expan- 
sion, which in turn ends at a transition to the expansion 
of Standard Big Bang cosmology We choose the initial 
conditions for the density of regular matter and for the 
value of (f> such that the temperature at which the Ekpy- 
rotic phase begins is higher than that at the time of equal 
matter and radiation in the Standard Big Bang expand- 
ing phase. In this way, we ensure that initial vacuum 
perturbations develop into a scale-invariant spectrum of 
cosmological fluctuations on all scales which are currently 
probed. 

In the following we study the evolution of the 
anisotropy in each of the periods of cosmological evo- 
lution mentioned above. 



A. Matter contraction 

We start by considering the period when the universe 
is dominated by a pressureless matter fluid, i.e., when 
the background equation of state parameter is roughly 
W = 0. In this phase, the mean scale factor evolves as 



a(t) 



2/3 



(26) 



where t E denotes the final moment of matter contraction 
and the beginning of the Ekpyrotic phase, and a E is the 
value of the mean scale factor at the time t E . In the 
above, t E is an integration constant which is introduced 
to match the mean Hubble parameter continuously at the 
time t„, i.e., 



Integrating the above yields the following expressions 
for the anisotropy factors 



= / ei{t')&t' = - 



3alH„ t-L 



2alM 8t 
3a 3 W 



-H(t). 
(30) 

From Eq. (30 1, we can observe that the anisotropy fac- 
tors 6i approach zero in the limit of t — > — oo. As the 
universe contracts with a matter-dominated equation of 
state, the absolute value of H increases and thus 9i be- 
comes larger as well. 

Before the Ekpyrotic phase, the universe is dominated 
by the pressureless matter fluid and the energy density 
of normal matter fluid evolves as p m ~ 3M 2 H 2 . How- 



ever, Eq. ( 29 1 shows that the effective energy density of 
anisotropics evolves as pg ~ H 4 and will therefore come 
to dominate. This is nothing but another way of stating 
the BKL instability in a contracting matter-dominated 
universe. In the following we will see that the Ekpyrotic 
phase of contraction cures this problem. However, in or- 
der to be able to enter this phase, the initial anisotropy 
cannot be too large for the model not to break down al- 
together. The requirement is that pg is less than p m at 
the transition time t , which implies 



(31) 



a requirement which we implement in our numerical dis- 
cussion below. 

Let us now move on to the discussion of the evolution 
of the anisotropy factors during the Ekpyrotic phase of 
contraction. 



B. Ekpyrotic contraction 



tw — ^TT 



3H W 



(27) 



Correspondingly, the mean Hubble parameter during the 
period of matter contraction is given by 



H(t) 



(28) 



Following Eq. ( 23 1 , one can immediately write down 



the time derivatives of the anisotropy factors as follows 



i(t)=M $i 



1 aim ' 



(29) 



where H E is the value of mean Hubble parameter at the 
moment t E . Note that the dimensional parameters of 
anisotropy M$ t i are therefore related to the values of 9i 
at the moment t E ; a smooth bounce will thus take place 
only if these values satisfy some constraints which we 
derive below. 



We assume a homogeneous scalar field </> which is ini- 
tially placed in the regime of <f> <C — 1 in the phase of 
matter contraction. In this case, the Lagrangian for <p ap- 
proaches the conventional canonical form and thus yields 
an attractor solution which is given by 



In 



2V r (t-f B _) 



(32) 



where t B _ is an integration constant which is chosen in 
order that the mean Hubble parameter at the end of the 
phase of Ekpyrotic contraction matches with the one at 
the beginning of the bouncing phase. This attractor so- 
lution corresponds to an effective equation of state 
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1 



2 

37 



(33) 



During the phase of Ekpyrotic contraction, the mean 
scale factor evolves as 



o(t) 



t-tr, 



(34) 
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where a B _ is the value of mean scale factor at the time t B _ 
which corresponds to the end of Ekpyrotic contraction 
and the beginning of the bouncing phase. Therefore, the 
mean Hubble parameter is given by 



H(t) 



<1 



t - U 



(35) 



where, in order to make H(t) continuous at the time i E 
one must set 



(36) 



Additionally, we require the mean scale factor to evolve 
smoothly and continuously at the time t E . This leads to 
the relation 



g B - 



(37) 



Given these preliminaries, we find that the time deriva- 
tives of the anisotropy factors evolve as 



i(t) 



Ma 



n3 



H 

Hi 



3q 



(38) 



and thus the effective energy density of anisotropy evolves 
as pg ~ H 6q . Hence, whereas pg still increases during the 
phase of Ekpyrotic contraction, the growth rate is much 
slower than in the matter-dominated phase (given the 
small value of q). In particular, the energy density in <f> 
increases much faster than that due to the anisotropics: 
the BKL instability is tamed in this manner by the Ekpy- 
rotic contraction phase. 

Integrating the above expressions for 8i we obtain 



»(*) 



i(t')dt' 



(2-3g)ffl3M fl ,, 
3(1 -3q)alH E 



-1 



3q 



2-3q 



H 



1-3?' 



(39) 



From the expression ( 39 1 we see that the absolute values 



of Pi see monotonically increasing when cosmic time t 
evolves from t E to the onset of the bouncing phase at time 
t B _ . However, they very rapidly converge to their limit- 
ing values, and since any constant part of the anisotropy 
factors can be absorbed in a rescaling of the spatial co- 
ordinates, it turns out that the actual anisotropy, effec- 
tively measured as the distance to an effective FLRW 
space, is effectively decreasing as the amplitude of the 
mean Hubble parameter increases: the Ekpyrotic phase 
thus provides a simple solution to the BKL instability 
growth in a contracting phase. 



C. Bounce phase 

In our model the scalar field evolves monotonically 
from cf) <C — 1 to <f> 3> 1. For values of cf> between 



/p/2 ln(2<7o) and 



/^721n(2 ff0 )/6 s 



ing one term in the denominator of g(<fi) dominates over 
the other at each transition time) , the value of the func- 
tion g((f>) becomes larger than unity and thus the universe 
enters a ghost condensate state. The occurrence of the 
ghost condensate naturally yields a short period of null 
energy condition violation and this in turn gives rise to 
a nonsingular bounce. 

As shown in Ref. |38J, we have two useful parameteri- 
zations to describe the evolution of the scale factor in the 
bounce phase. One is the linear parametrization of the 
mean Hubble parameter 

H(t) ~ Tt, (40) 

and the other is the evolution of the background scalar 



-t 2 /T 2 



(41) 



where the coefficient T is set by the detailed microphysics 
of the bounce. The coefficient T can be determined by 
matching the detailed evolution of the scalar field at the 
beginning or the end of the bouncing phase, which will be 
addressed in next subsection. Thus, during the bounce 
the mean scale factor evolves as 



a(t) 



a = e 2 



t/- 



(42) 



Note that a nonsingular bounce requires that the total 
energy density vanishes at the bounce point. The total 
energy density includes the contributions from the mat- 
ter fields and the anisotropy factors. This leads to the 
following result for the value of <f) B 



(go - mi 

3/3 

2(ffo - 1) 

3/3 




12P(V + Pm + pg) 



(.90 - l) 2 M p 4 



M 



(43) 



where we have made use of approximations that p m and 
pg are much less than Vq and Vq <C Ml in the second 
line. These approximations must be valid for the model 
to hold since both p m and pg are greatly diluted in Ekpy- 
rotic phase and Vq is the maximal absolute value of the 
potential of (f> which, according to the observational con- 
straint from the amplitude of cosmological perturbations, 
must be far below the Planck scale. 

Now we calculate the anisotropy factors in the bounc- 
ing phase. We first study their time derivatives, which 
are given by 



,(t) ~ Mg^e- 



(44) 



Thus we see that the integration constants Mgj can 
be interpreted as the values of 9i at the bounce point 
t B =0, and they are the maximal values Qi will ever 
take throughout the whole cosmic evolution. It is again 
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easy to perform the integrals and obtain the anisotropy 
factors: 



phase we find the following approximate solution for the 
evolution of d>: 



i(t) 



. (45) 



where Erf is the Gauss error function. Around the 



bounce point, we can perform a Taylor expansion of (45) 



up to leading order and thus obtain the following approx- 
imate expression, 



i(t) 



+ M e<i {t-t B _) 



(46) 



which is a linear function of cosmic time. We find that 
the growth of the anisotropy at linear order only depends 
on the coefficients Mg^. However, at nonlinear order in 
cosmic time, the growth would also depend on the slope 
of the bounce phase characterized by the parameter T. 
Finally, we have the relations 



a B e- 



(47) 
(48) 



by matching the mean scale factor and the mean Hubble 
parameter at the beginning of the bouncing phase. 



D. Fast-roll expanding phase 

After the bounce, the universe enters the expanding 
phase, which typically begins with a period of fast-roll 
expansion for the background during which the universe 
is still dominated by the scalar field <f>. During this stage, 
the motion of (f> is dominated by its kinetic term while the 
potential is negligible. Thus, the background equation of 
state parameter is w ~ 1. Correspondingly, the mean 
scale factor evolves as 



a(t) ~ a B+ ( V 



1/3 



(49) 



<t>(t) - & 



+ a 3 (i) 



,e-^'fa, (52) 



where we have applied (41 1 in the second equality. This 
implies that 



< 12 



M 2 J 2 B H 2 
2e 2t B + / T2 H* 



(53) 



Moreover, the Friedmann equation requires that p^ 
3M 2 H 2 in the fast-roll phase, so that T 2 is given by 



T 2 



2H 2 



T 2 In 



9/3# 2 + 



(54) 



During the fast-roll expansion epoch, one can solve for 
the time derivatives of the anisotropy factors 



i(t) =s M e , 



(55) 



from which it follows that the effective energy density 
of anisotropy evolves as pg — pg yB+ H 2 / H 2 + , and thus 
its ratio relative to the total energy density does not 
change (recall that the absolute value of this ratio is also 
very small because of the decrease of the shear contri- 
bution during the Ekpyrotic phase of contraction dis- 
cussed above). On the other hand, the matter fluid, 
having u> m = 0, sees its energy density going as p m — 
p m , B+ H/ H B+ . and therefore its contribution to the back- 
ground will catch up with that of the scalar field and the 
anisotropy, bringing the fast roll phase to an end and ini- 
tiating the transition to the usual expansion history of 
Standard Big Bang cosmology. 

The anisotropy factors in the fast roll phase are given 

by 



where t B+ represents the end of the bounce phase and 
the beginning of the fast-roll period, and a B+ is the value 
of the mean scale factor at that moment. Then one can 
write down the mean Hubble parameter in the fast-roll 
phase 



H(t) 



1 



3(t-t B+ ) 



(50) 



and the continuity of the mean Hubble parameter at t B , 
yields 



- _ 1 



(51) 



Recall that, in Eq. (41), we made use of a Gaussian 



parametrization of the scalar field evolution in the bounce 
phase, with characteristic timescale T. In the fast roll 



i(t) 



alM e ,i 



In 



t-t 



B + 



(56) 



so they are safely growing only logarithmically during 
this phase. 

To conclude this section, it is clear that the anisotropy 
contribution can easily be made to remain small through- 
out the entire evolution, i.e. including the matter con- 
traction and the bounce phases. We now turn to the 
actual constraints that should be imposed on the param- 
eters of our model. 



E. Theoretical constraints 

In this subsection, we study the theoretical constraints 
on the scenario of anisotropic nonsingular bounce. It is 
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convenient to introduce an effective e-folding number of 
the Ekpyrotic phase as follows 



AC = In 



= (1-9) In 



(57) 



which characterizes the variation of the mean value of 
the conformal Hubble scale aH during this period. 

Recall that the anisotropy is not allowed to dominate 
before the onset of the Ekpyrotic phase. This is the con- 
Applying (37 1 and (57 1, 



straint 
derive t 



(311 



we can further 



re following relation for the anisotropy at the 
bounce time: 



E 



Mi 



-2(l-3g)jV E /(l-g) 



(58) 



which is. as expected, exponentially suppressed by the 
number of e-foldings of Ekpyrotic contraction. Note that 
one can introduce a relative density parameter to de- 
scribe the contribution of the anisotropics as follows 



ZMIH 2 ' 



(59) 



This parameter must be less than the amplitude of the 
observed anisotropics in the cosmic microwave back- 
ground (CMB), which is of order f2 port ~ 0(1O~ 5 ), 
hereby defining f2 per t as the amount of energy density in 
perturbations, as given by CMB observations. This leads 
to the following constraint on the number of e-foldings 



1^ Q H 2 



< n 



pert i 



(60) 



where we have made use of ( 55 ) . Provided that this 



constraint is satisfied, the current anisotropy will be be- 
low the current upper bound provided that it was ini- 
tially small enough not to prevent the onset of Ekpyrotic 
contraction. Combining this inequality and the approxi- 
mate relation ( 58 I , we can obtain a lower bound on the 



e-folding number of the Ekpyrotic phase 



AC > 



i-q 

2(1 -3q) 



In 



1 



pert 



(61) 



Note that the nonsingular bounce model studied in the 
present paper belongs to "fast bounce" category and the 
absolute values of H B and H B are of the same order. 
Thus, we can simply consider H B to be the energy scale 
of the bounce. As an example, if we choose q — 0.1 and 
^pert = 10~ 5 , then the relation ([61]) yields AC > 7.4. 
This result implies that the Ekpyrotic phase is rather 
efficient at lowering the contribution of anisotropics since 
a mere few e-folds of Ekpyrotic contraction are enough 
to damp any reasonable initial anisotropy to negligible 
level. 

Note that the e-folding number AC can also be con- 
strained on observable scales of CMB experiments, and 



we expect this may impose a much more stringent con- 
straint on the anisotropy parameters. For that, one 
could require that the anisotropy contribution be smaller 
that the observed level of perturbations, namely f2 per t ~ 

io- 10 . 

It is interesting to note that increasing the Ekpyrotic e- 
fold number to, say, AC ~ 15 while keeping q ~ 0.1, leads 

10 10 



to relation (31) to yield fig 



O 



pert 



pert ■ 



This means that after a sufficiently long Ekpyrotic con- 
traction, the anisotropy contribution is totally negligible, 
even compared with the amplitude of primordial pertur- 
bations. At this particular level however, the anisotropy 
contribution could be comparable to second order in 
terms of the primordial curvature perturbation expan- 
sion, and thus could contribute to the primordial non- 
Gaussianities of local shape with / NL ~ 0(1). 

If one wants to make successful contact with late time 
cosmology, there is a second constraint that should be im- 
plemented on the model, namely that the fast roll phase 
ends before the time of Big Bang Nucleosynthesis (BBN). 
Roughly speaking, the energy density of regular matter 
does not change much during the phase of Ekpyrotic con- 
traction (for small values of q). On the other hand, the 
density of <p grows rapidly. In the fast roll phase of ex- 
pansion, the decrease in the density of regular matter is 
no longer negligible. Hence, the energy density of matter 
will be much lower at the time t F when p m (t F ) — p<p{t F ) 
than at the time t F when Ekpyrotic contraction begins. 
In fact, it is straightforward to derive that 



H„ < \H„ 



-(l-3 9 )JV E /(l- 9 ) 



(62) 



showing that the value of H F should be much less than 
\H E \. The Hubble rate H F is associated with the ini- 
tial temperature T F when the expansion begins to follow 
the Standard Big Bang evolution (it is the equivalent of 
the temperature of reheating in inflationary cosmology). 
Specifically, the relation is 



1/2 



7TT 2 



9.5M D 



(63) 



where g s is the effective partible number for radiation. 
As a consequence, in analogy with inflationary cosmol- 
ogy, the constraint ( 62 1 leads to an upper bound on the 



effective "reheating" temperature: 



T < 



3M P1 |g B 

1/2 



-(2-3 9 )Af B /[2(l- g )] 



(64) 



in our nonsingular bounce model. From the BBN con- 
straint, we find that the lower limit of the "reheating" 
temperature is of the order 0(MeV). If we consider this 
lower bound and take g s ~ 100, AC ~ 30 and q ~ 0.1, 
then we find that H B+ > 10~ 17 M P1 which can easily be 
implemented in the model, as we shall see in the following 
numerical calculations. 
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Numerical estimates 



To illustrate that a nonsingular bounce can be achieved 
in our model, we numerically solved the background 
equations of motion. Expressing all relevant functions 
and parameters in the corresponding units of the reduced 
Planck mass M pl , we set 

V = 10- 7 , 9o = l.l, /3 = 5, 7=10~ 3 , 
b v = 5, b g = 0.5, p = 0.01, q = 0.1 (65) 

to illustrate the calculations. 

Moreover, we consider the following parameters of the 
matter fluid and the anisotropy 



= 2.8 x 1CT 10 , M e ,i = 2.2 x 1CT 6 , 
= 3.4 x 1CT 6 , M e 3 = -5.6 x 1CT 6 , 



(66) 



and choose as the initial conditions for the scalar field 
the following: 



= -2, 



= 7.8 x 10" 



(67) 



The actual computation also requires the initial value 
of the mean Hubble parameter, which is determined by 
imposing the Hamiltonian constraint equation. Figs. [2] 
and [3] show the evolution of the Hubble parameters and 
"effective" energy densities for matter components and 
for the anisotropy, respectively. 

From Fig. [2j one can see that Hubble parameters 
along all spatial coordinates evolve smoothly through the 
bouncing point with an approximate dependence on cos- 
mic time which is linear. The maximal value of the mean 
Hubble parameter, which we denote as the bounce scale 
H B , is mainly determined by the value of the potential 
parameter Vq- Specifically, H u is of order 0(1O~ 4 M PI ) in 
our numerical result. We also note that the bounces oc- 
curring in the three spatial directions do not occur at ex- 
actly the same moment - a consequence of the existence 
of anisotropy. This could leave a smoking gun signature 
for detecting nonsingular bounce cosmology in high accu- 
racy CMB experiments since the difference in the times 
of the bounces along various spatial coordinates would 
affect the ultraviolet (UV) modes of primordial pertur- 
bations passing through the bouncing phase. We leave 
this issue for a forthcoming investigation. 

From Fig. [3] one can easily see that the universe in 
our model experiences four phases, which are matter con- 
traction, Ekpyrotic contraction, the bounce, and fast-roll 
expansion, in turn. At the beginning, the universe is 
dominated by the matter fluid. At some point (time t E ) 
during the phase of contraction, the contribution of the 
scalar field becomes dominant, and the universe enters 
the Ekpyrotic phase. Note that the effective energy den- 
sity of anisotropics grows faster than p m but slower than 
Ptf, during the matter contraction phase. Thus, if p@i 
does not dominate over the background before i E , it will 
never become dominant throughout the whole evolution, 
as already discussed in the previous sections. After the 
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0.0 
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0.0 2.0x1 4 
t 
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Figure 2: Time evolution of the Hubble parameters H (black 
line) and Hi (red dashed, blue dotted and magenta dot-dashed 
lines for the Hubble expansion rates along the x%, X2, and X3 
axes, respectively), in units of the reduced Planck mass M p , , 
with background parameters given by Eqs. (65 1 and (66 1, 
and initial conditions as in (67 1. The main plot shows that 



a nonsingular bounce occurs, and that the time scale of the 
bounce is short (it is a "fast bounce" model). The inner insert 
shows a blowup of the smooth Hubble parameters during the 
bounce phase: this zoomed-in view of the Hubble parameters 
around the bounce point shows that the Hubble rates vanish 
at different times, so that the scale factors bounce at different 
times as well. 




-3.0x1 4 0.0 3.0x1 4 6.0x1 4 9.0x1 4 
t 



Figure 3: Time evolution of the "Effective" energy densities of 
the scalar field (full black line), the matter fluid p m (dot- 
dashed green line) and the anisotropy factors pei (red dashed, 
blue dotted and magenta dot-dashed lines), with same initial 
conditions and background parameters are in Fig. [2] 



bounce, the scalar field <j) enters a fast roll phase with 
an effective equation of state equal to unity. As a con- 
sequence, the energy densities p^ and pei dilute at the 
same rate, and finally the matter fluid catches up with 
the density of </> at the time t F . 

Fig. [4] shows the evolution of the anisotropy factors 
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Figure 4: Time evolution of the anisotropy factors 9i (lower 
panel) and their time derivatives 9i (upper panel), with same 
initial conditions and background parameters are in Fig. [2] 
The anisotropics increase during the contracting phase but 
rapidly approach constant values in the following expanding 
phase. 

9i and their time derivatives. Although the anisotropy 
functions grow during the contraction, they evolve to- 
wards constant values in the expanding epoch. There- 
fore, after the time t F , one can rescale all scale factors by 
absorbing the asymptotic factors in a redefinition of the 
coordinates, and we finally get an isotropic universe. It 
implies that at the level of homogeneous cosmology the 
anisotropies do not destabilize our nonsingular bounce 
model. This can also be read from the upper panel of 
Fig. [i] which shows that 9\ approach zero after a suffi- 
ciently long period of expansion. 

In order to better characterize the anisotropy quanti- 
tatively, we can define the so-called shear parameters 

a, = ^e 2e % (68) 

and the density parameters 

n t = =^-, (69) 

where the subscript "1" represents </>, m and 9, respec- 
tively. Fig. [5] shows the numerical solution we obtained 
for their time development. The shear functions increase 
up to their maximal values at the bounce point, after 
which they rapidly decrease to end up vanishingly small 
when the universe connects with the Standard Big Bang 
evolution. From the evolution of the density parameters, 
we see that the contribution of the anisotropy only grows 
relative to the dominant density in the phase of matter 
contraction, but it then rapidly decreases in the Ekpy- 
rotic phase and in the fast roll phase. 

Note that all numerical calculations shown here are 
meant to illustrate the discussion of the previous sections. 



Figure 5: Time evolution of the density parameters fl^, f2 m , 
and 0,0 (upper panel), and of the shear function <7; (lower 
panel), with same initial conditions and background parame- 
ters are in Fig. [2] 

Indeed, the parameters chosen do not satisfy the bounds 
imposed by CMB observations, so the effects of setting 
non vanishing initial anisotropies can be enlarged so as 
to be visible at all. Assuming parameter values taking 
into account the experimental constraints would lead to 
figures exactly similar to those obtained for a regular 
isotropic bouncing model. 

IV. CONCLUSION AND DISCUSSION 

We have studied homogeneous but anisotropic cos- 
mological solutions of the Galileon type action with an 
Ekpyrotic scalar field potential, as originally introduced 
in [38J . The model was constructed to yield non-singular 
isotropic bouncing cosmological solutions. In general, the 
contracting phase of a bouncing cosmology suffers from 
a BKL instability against the growth of anisotropies. By 
studying the solutions of Bianchi-type homogeneous so- 
lutions in which the scale factors in different spatial di- 
mensions are allowed to be independent we have stud- 
ied the evolution of the anisotropies in our model. We 
have shown that the anisotropies are damped in the phase 
of Ekpyrotic contraction, as was already conjectured on 
qualitative grounds in |38| . They remain small during the 
non-singular bounce phase, and are further suppressed 
once space begins to expand again. This (counter) ex- 
ample shows that a no-go theorem concerning smooth 
bounces cannot be formulated, at least in terms of BKL 
instabilities. 

It would be of interest to study the evolution of lin- 
earized inhomogeneities in our anisotropic background 
and to see if there are any late time signatures of 
the initial (pre-Ekpyrotic) contracting phase when the 
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anisotropics might be quite large. This is left for future 
work. 
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